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Abstract 

Delsarte's method and its extensions allow to consider the upper bound 
problem for codes in 2-point-homogeneous spaces as a linear programming 
problem with perhaps infinitely many variables, which are the distance 
distribution. We show that using as variables power sums of distances this 
problem can be considered as a finite semidefinite programming problem. 
This method allows to improve some linear programming upper bounds. 
In particular we obtain new bounds of one-sided kissing numbers. 



1 Introduction 

The Delsarte linear programming method is widely used for finding bounds for 
error-correcting codes, constant weight codes, spherical codes, sphere packings 
and other packing problems in 2-point-homogeneous spaces [H [H [TOl [121 IE] ■ 
Linear Programming (LP) is a special case of Semidefinite Programming (SDP) 
which deals with optimization problems over symmetric positive semidefinite 
matrix variables with linear cost function and linear constraints. 

Recently, Schrijver [SJ using SDP improved some upper bounds on binary 
codes. Even more recently, Schrijver's method has been adapted for non-binary 
codes (Gijswijt, Schrijver, and Tanaka [TT]) and for spherical codes (Bachoc 
and Vallentin [5]). (See also our paper In fact, this method using the 

stabilizer subgroup of the isometry group derives new positive semidefinite con- 
straints which are stronger than the linear inequalities in the Delsarte linear 
programming method. 

In this paper we define semidefinite programs whose optimal solutions give 
upper bounds for codes in a 2-point-homogenous space and strengthen the LP 
bounds. Our approach is based on two ideas: to consider power sums of dis- 
tances instead of the distance distribution, and to use the positive-semidefinite 
property of zonal spherical functions. In fact, this method is different from 
Schrijver's method and can be applied for any 2-point-homogenous space. 

The paper is organized as follows: Section 2 reviews on the LP method 
for codes in 2-point-homogenous spaces. Section 3 recalls a theorem due to 
Sylvester. Section 4 introduces the SDPq bounds for codes and shows that 
these bounds are equivalent to the LP bounds. Section 5 considers examples of 
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applications of SDPq . Section 6 shows that some recent extensions of Delsarte's 
method can be reformulated as SDP problems (SDPA). Section 7 extends the 
SDPA bounds to subsets of a 2-point-homogeneous space and shows that some 
upper bounds for codes can be improved. 

Our computations were performed by the program SDPA-M, version 2.0 (see 
|http: / /sdpa.is.titech.ac.jp/SDPA-M / ) . After solving the SDP with SDPA-M 
we checked independently whether the solution satisfies the desired constraints. 
This can be done using rational arithmetic only. 

2 The linear programming bounds for codes 

The main goal of the present section is to give a brief introduction to zonal 
spherical functions and the linear programming bounds for codes. We are based 
on the simplified account, given by Sloane [51 Chapter 9], of the general ma- 
chinery developed by Kabatiansky and Levenshtein [1^ . Here we are using the 
same notations as in [5] and we hope that the reader will have no difficulty using 
Sloane's exposition to fill up details and proofs. 

We say that a G-space M is a 2-point-homogeneous space if (i) M is a metric 
space with a distance p defined on it; (ii) M is a set on which a group G acts; (iii) 
p is strongly invariant under G, i.e. for x,x' ,y,y' G M with p{x,y) = p{x' ,y') 
there is an element g G G such that g{x) = x' and g{y) — y' . 

These assumptions are quite restrictive. In fact, if G is infinite and M is a 
compact space, then Wang [30] has proved that M is a sphere; real, complex 
or quaternionic projective space; or the Cayley projective plane. However, the 
finite 2-point-homogeneous spaces have note yet been completely classified (see 
for the most important examples and references [S]). 

With any 2-point-homogeneous space M and an integer number A: > are 
associated the zonal spherical function $fc(i) such that {$fc(t)}fc=o,i,2,... are or- 
thogonal on T := {t(x,?/) : x,y M}, where r is the certain function in p 
( i.e T{x,y) — F{p{x,y))) defined by M. For all continuous compact M and 
for all currently known finite cases: ^k{t) is a polynomial of degree k. The 
normalization is given by the rule: ^k{To) — Ij where tq := t{x,x). Then 

$o(0 = 1. 

Let us consider two examples. 

Example 1: M = Hamming space with t{x, y) — p{x, j/)=Hamming distance. 
Then tq — 0. Here $fc(i) is the Krawtchouk polynomial Kk{t, n): 

*'<"-^'*"'K:)"S<-K:)(:::)- 

Example 2 (the most extensively studied case): M = unit sphere S"~^ C R" 
with t(x, y) = cos p(x, y) = (x, y), where p{x, y) is the angular distance between 
X and y. In this case tq = 1. The corresponding zonal spherical function ^^(i) 
is the Gegenbauer (or ultraspherical) polynomial G^^^t). 
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There are many ways to define the Gegenbauer polynomial G^^^ (see [H 
[TOl [T2I [T8l [23l [26]). which is a special case of Jacobi polynomials Pjf'^^ with 
a = /3 = (n — 3)/2 and normalization G^"^(l) = 1. G^"^ are orthogonal on 
[—1, 1] with weight (1 — t2-j(n-3)/2_ ^jgQ Qi^) (.g^j^ ]-,g defined by the recurrence 
formula: 

Mn) (2fc + » - 4) ^ - (fc - 1) Gi% 
Go -1, Gi -t, ...,G, - k + n-3 " 

The main property for zonal spherical functions is called "positive-definite 
degenerate kernels" or "p.d.k" [5]. This property first was discovered by Bochner 
[6] (general spaces) and independently for spherical spaces by Schoenberg [26] . 

Now we explain what the p.d.k. property means for finite subsets in M. 

Theorem 2.1 ((nHHHHl). Let M &e a 2-point-homogeneous space. Then for 
any integer k > and for any finite C — {xi} C M the matrix ($fc(T(xi, Xj))) 
is positive semidefinite. 

This theorem implies the fact that plays a crucial role for the linear pro- 
gramming bounds. For any positive semidefinite matrix the sum of its entries 
is nonnegative. Then 

Theorem 2.2 ( [101 [121 HI] ) • For any finite C = {xj C M we have 

|c| \c\ 

^^<i>fc(r(x„a;j))>0. 

i=i j=i 

Since ^kiT^Xi, Xi)) ~ $fe(To) = 1, this theorem implies 
1 

^1 



^kiT{x,,Xj)) > -1, A; = 0, 1, 2, ... . (1) 



Let 5* be a fixed subset of T. We say that a finite subset G C M is an S-code 
if T{x,y) G S for all x,y & C,x ^ y. The largest cardinality |G| of an S'-code 
will be denoted by A{M, S). 

Usually, 5* — T n [a,b] for some interval [a,b]. For instance, if we take 
S = Tr\[d,n] for Example 1, then an S'-code is an error-correcting code of 
minimal distance d. For Example 2, if S" = [— l,cos6'], an 5-code (or a 6'-code) 
is a spherical code with minimal angular distance 9. 

The distance distribution {a*} of G is defined by 

at := T-^ (number of ordered pairs x,y d C with t{x, y) =t). 
|G| 

We obviously have 

aro = 1, at > 0, t e T, ^ at = |G|. (2) 

teT 
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(1) and (2) make it possible to regard the problem of bounding A(M, S) as 
a linear programming problem: 

Primal problem (LPP): Choose a natm'al number s, a subset {ti, . . . , r,} of S, 
and real numbers , • ■ • , otr^ so as to 

maximize a,-! + . . . + 

subject to 

s 

"r, > 0, i = l,...,S, ^ar,$fc(Ti) > -1, = 0, 1, . . . . 
i=l 

This is a linear programming problem with perhaps infinitely many un- 
knowns at and constraints (f ), (2). If C is an 5-code then its distance distribu- 
tion certainly satisfied the constraints (1), (2). So if the maximal value of the 
sum a^-i + • • • + ttr, that can be attained is A* , then A(M, S) < 1 + A* . (The 
extra 1 arises because the term am = 1 doesn't occur in this sum.) 

Dual problem (LPD): Choose a natural number N and real numbers /i, . . . , /jv 
so as to 

minimize /i + . . . + /^r 

subject to 

N 

h>Q, k^i,...,N, 5]/fe$feW<-i, teS. 
fc=i 

Thus, we have 

Theorem 2.3. // A* is the optimal solution to either of the primal or dual 
problems, then A(M, S)<1 + A*. 

Note that for the dual problem any feasible solution is an upper bound to the 
optimal solution of the primal problem. Using this (or using directly Theorem 
2.2) it is not hard to prove: 

Theorem 2.4 ([SJ [TOl [HI HH [2I])- Suppose that f is an nonnegative linear 
combination of {^k}k=o,i-- -,N for some N, i.e. 

fit) = foMt) + fMt) + ■■■ + /iv$iv(i), where fk > 0. 
If fit) < for all teS and fo > 0, then 

AiM,S) < 

fo 

This theorem and its extensions have numerous applications in coding theory 
and sphere packings (see [5] for references, and for some recent achievements 
[15l[23l[T6l[T7l[T8l|T9]). Since one of the most important steps in this approach 
was taken by Delsarte [9] it is often called Delsarte's linear programming bound 
or Delsarte's method. 
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3 Sylvester's Theorem 



One of the classical problems for polynomials is the problem of finding the 
number of positive roots of a real polynomial. There are several fundamental 
facts in this field (see [231 Chapter 1]). One of the powerful methods is called 
Hermite's method. This method is based on (Hankel's) quadratic forms which 
coefficients depend on power sums Sk := + . . . + tj^, where ti are roots of a 
polynomial. Actually, this method gives an answer for the following question: 
Suppose we know for complex numbers ii, . . . , i„ only its power sums s^, where 
Sk are real numbers. How to determine the number of real and positive (or 
nonnegative) ti? 

Theorem 3.1 (Sylvester, [H Section 1.4]). 

(i) Suppose that for complex numbers ii, . . . , t„ all power sums Sk are real num- 
bers. Then the number of distinct real ti is equal to the signature of the quadratic 
form defined by the matrix 



( 



R„ 



So 

Sl 



Si 
S2 



\ Sn~l 



Sn-1 \ 
Sn 

S2n-2 / 



(ii) If all ti are real, then all of them are nonnegative if and only if the matrix 
Fn is positive semidefinite, where 



Fr,. 



/ s. 



S2 
S3 

Sn+1 



Sn \ 

Sn+1 



S2r, 



J 



So we have that all ti are real and nonnegative if and only if i?„ and Fn are 
positive semidefinite matrices (i?„ >: 0, Fn >: 0), and the number of distinct 
ti equals rank(i?„). It is clear that Sylvester's theorem gives an answer for the 
same question when all ti are real and belong to some interval [a,b]. 

Now we introduce the following matrices: 



I Sl - aso 
S2 — asi 



52 — as\ 

53 - aS2 



\ Sni aSni—1 Sni-\-l OSn 



F^nih) 



( bso - Sl 

bs\ — S2 



bsi — S2 
bs2 - S3 



\ bSni—l Sjn bSjn S^n+l 



Sm ^Sni^l \ 
Sm+1 ^Sni 

S2m-1 — aS2m-2 / 

bSni—1 Sni \ 
bs^n Sni-\~1 

bS2m-2 — S2m-1 / 



5 



and 

/ Rm 

Fm(a, &) =iJm(so,si...,S2m-i,[a,6]) := F+{a) 

V F-{b) 

Theorem 3.2. Consider real numbers ii,...,i„ in [a.b]. Then for any nat- 
ural number m the matrix H^(a,b) is positive semidefinite. Moreover, for 
m > rank(i?„) the converse holds, i.e. if Hm(a,b) >: 0, then complex num- 
bers ti,...,tn with real are real numbers in [a, b] . 

Proof. The proof is based on the classical Hermite's proof of Sylvester's theorem. 
Let 

n 

^2 



G(U, V, W) = ^ (Ui + tkU2 + ...+ f^-^Umf 

fc=l 

n n 
■Y,{ik-a){v^ + ... + t^-\m? + ^ (fo - + ■ • • + tT^Wmf 

fe=l fe=l 



^,3 1,3 

It is easy to see that pij = Si+j-2, Qij = Si+j-i — asi+j-2, rij = bsi 
Si+j_i. Therefore, H^{a,b) is the matrix of the quadratic form G. If all ti G 
[a, b] , then G is nonnegative for any u, v, w, i.e. G ^ 0. 

Note that r = rank(i?„) is the number of distinct ti = the total number of 
squares of linear forms (positive and negative) in the quadratic form G(u, 0, 0). 
Then G(u, 0, 0) )^ yields that all ti are real numbers. Moreover, G(0, v, 0) ^ 
and G(0, 0, w) ^0 imply that all of them lie in [a, b]. □ 

In fact, the constraint Hm{a, b) ^ doesn't depend on n = sq. Indeed, let 

:= 1±111±^ = £fc, H^{a,b) = -Hm{a,b). 
n n n 

In other words, Hm{a, b) can be obtained by substituting Sk for Sk in Hm{a, b): 

Hm{a, b) = -ffm(si, . . . , S2TO-1, [a, b]) := -ffm(l, Si,..., S2m-1, [a, b]). 

Thus Hm{a, 6) t if and only if Hmia, b) h 0. 

Let us define in (2m — 1) - dimensional linear space {(si, . . . , S2m-i)} — 
j^2)n-i domain 

A„(a,6) := {{xi, . . . , X2,n-i) e R^™"^ : H^{xu . . . , X2m-i, [a,b]) h 0}. 

Theorem 3.3. For any real numbers a,b, a < b, Am,{a, b) is a compact convex 
set in R^"*"^ of dimensions 2m — 1. 

Proof. Since — oo < a < b < +oo and \xk\ < maxdaj*^, we see that 
Am{a,b) is a compact set. Note that Hm depends on xi, . . . ,X2m-i linearly. 
Clearly, the convexity holds for any positive semidefinite matrix with entries 
that are linearly dependent on the variables xi, . . . , X2m-i- Since Hm{a, 6) ^ 
define a nonempty open set, we have dim(Am(a, b)) = 2m — 1. □ 
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4 LP and SDPq bounds 

In this section we consider the hnear programming bound {LP) and the simplest 
version of the SDP bounds (we denote it by SDPq) for codes. 

From here on we assume that is a polynomial of degree k, $fe(ro) — 1, 

and S — T n [a,b] (the most interesting case for coding theory and sphere 
packings). 

In fact, the optimal solution A* of the LPP and LPD problems in Section 
2 depends only on the family of polynomials $ := {^k{t)}k=o,i,... and [a, b]. We 
denote 1 + ^* by LP{<P, [a,b]). 

Since ^k{t) is a polynomial of degree k, we have: 

$fc(i) = Pko + Pkit + ...+ Pkkt^- 

Then 

k 

$fc(ii) + . . . + ^k{tt) = Y^PkdSd, sd = ti + ...+t: 



d 



i=0 



Let C = {vi} be an S'-code on M, and let Tij = T{vi,Vj). Note that the 
number of ordered pairs {vi,Vj), i j, equals n = |C|(|C| — 1). Then (1) can 
be written in the form: 

k 

y + PkO + X! PkdXd > 0, (3) 



d=l 



where 



1 - Sd rf 



From Theorem 3.2 we have for any m: 

H^{xi,...,X2,n^i,[a,b])yO. (4) 

A positive semidefinite constraints (3) and (4) make it possible to regard 
the problem of bounding ^(M, S) as a semidefinite programming problem. The 
standard form of the SDP problem is the following [551 US] ■ 

Primal Problem: 

minimize cixi + . . . + cex^ 

subject to 

X>zO, where X = Tixi + ... + Tixi - Tq. 

Dual Problem: 

maximize (Tq , Y) 

subject to 

(T„y) = c„ i = i,...i, 
y ^ 0. 
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Here Tq,Ti, . . . ,T£, X, and Y arc real Nx N symmetric matrices, (ci, . . . , q) 
is a cost vector, (xi, . . . ,Xi) is a variable vector, and by {A, B) we denote the 
inner product, i.e. {A,B) = Tr{AB) = '^aijbij. 

Now we introduce the SDPo{^, [a,b],m) bound. 

SDPo Problem: Choose a natural number m and real numbers y,xi, . . . , X2m-i 
so as to 

minimize y 

subject to 

k 

V + ^PkiXi > -Pko, k = l,...,2m-l, 

Hm{xi,. . .,X2m-i, [a,b]) h 0. 

Note that in (3): \C\ = (1 + y)/y. Thus 
Theorem 4.1. If y* is the optimal solution of the SDPq problem, then 

A{M,S) < SDPoi^,[a,b],m) := 1±^. 

y* 

Since we just substituted Hm >z for t G S in the LPP problem, we can 
expect that SDPo{^,[a,b]) = LP{^,[a,b]). It is not hard to give a formal 
proof that LPP = SDPq and LPD = Dual SDPq. We think that it becomes 
more clear by geometrical arguments. Actually, we are going to prove that 
SDPq = LPD. 

Consider in R^, N = 2m — 1, the domain defined by (3): 

k 

^m{y) ■■= {x e : "^PkiXi > -y-pko, k = l,.. .,N}. 

i=l 

Clearly, Am{y) is a convex cone in R^. 
Let 

k 

Uk = Pko + ^PkiXi, k=l,...,N. 

It is easy to see that (ui, . . . , ujv) is a nondegenerate coordinate system of R^. 
So in these coordinates 

Am.{y) = {{ui, ■ ■ ■,un) e R-'^ ■.Ui>-y for all i}, 

i.e. Am(y) is a pointed cone with vertex (— y, . . . , —y). 

For a given m the SDPq problem is asking to find the minimal y = y*{m) 
such that the intersection of Am{y) and Am{a, b) is not empty. Then y < y* if 
and only if the intersection of these convex sets is empty. 
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For < y < y* consider a separating hyperplane 

i(u) = /o + flUi + . . . + fNUN 

such that 

L(u) < 0, Vu e A„.(a, 6), L{u) > 0, Vu e A,„(y). 

Clearly, L(u) > yields fk > 0, fo > y(/i + - ■ ■ + fN)- On the other hand, it can 
be proved that L(u) < if and only if f {t) = fo + + . . . + /Ar"I>jv(<) < 

for all t e [a, b]. Then Theorem 2.4 yields I + A* = {I + y*)/y*. 
Thus we have proved: 

Theorem 4.2. S'L'Pol*, [a, &]) = i-P($, [a, b]). 

In fact, for a continuous M the LPP and LPD problems are not finite linear 
programming problems. These problems can be solved only via discretization. 
For instance, Odlyzko and Sloane [21] (= [8l Chapter 13]) apphed LPD for 
upper bounds on kissing numbers, where they replaced S by 1001 equidistant 
points in S. For the LPP problem it is not clear how to do a discretization 
of {ar}- On the other hand, for a given m the SDPq is a finite primal SDP 
problem. As a by-product of solutions of this problem we have bounds on \C\ 
and power sums Sk (see the next section). 

Remark. One of the important problems for SDPq is the problem of detecting 
an optimal m. We see that it is the same problem as the problem of optimal 
degree A'^ of a polynomial in Theorem 2.4. Our arguments show that m = 
r(Ar+l)/2l. 

The problem of finding optimal polynomials for Delsarte's linear program- 
ming bounds has been the subject of many papers. The most fundamental 
results in this problem were established by Levenshtein (see [15]). 

5 Bounds via SDPq 

In this section we consider three examples of spherical codes and define the main 
steps of SDPq. 

1. It was shown by Levenshtein |14j and independently by Odlyzko and Sloane 
|21| that the maximal number of nonoverlapping unit spheres in (resp. R^**) 
that can touch another unit sphere is 240 (resp. 196560). In other words, they 
proved that the maximal size of a spherical 7r/3-code in dimension 8 (resp. 24) 
is 240 (resp. 196560). Bannai and Sloane [4](=[8l Chapter 14]) proved that the 
maximal arrangements in these dimensions are unique up to isometry. 

Consider an application of the SDP method for these spherical codes. Let 
us suppose that we have no idea about the kissing number in eight dimensions 
and how a maximal kissing arrangement looks like. Now we try to find it by 
SDPo. 

Step 1: Semidefinite Programming. 
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Let us apply the SDPq with [a,b] = [-1,1/2] and $fc(i) = G[^>{t). For 
TO = 3,4,5 we obtain that c*(3) = 324, c*(4) = c*(5) = 240. (Here c*(to) = 
{l + y*{m))/y*{m) [see Section 4].) In fact, for any m > 4 we have c*(to) = 240. 
So we can expect that the maximal size of a kissing arrangement C — {vi} is 
240 and for to > 4: rank(i?m) = 4 (the number of distinct t,y = T{vi,Vj), see 
Section 3). The SDPq for to = 4 gives: 

xl=x;=xl^x*j^ -1/239, xl = 29/239, xl = 8/239, Xq = 11/956. 

Thus if there exists C with c = |C| = 240 and rank(i?,„) = 4, to > 4 , then 

4 

^ a, = ^fe, fc = 0, . . . , 7, where /3o = 239, := 239x1, k > I. 

1=1 

Step 2 : Feasible solution for the distance distribution. 

Let us consider the following problem: 
Choose four real numbers ti, — 1 < ti < t2 < is < ^4 < 1/2 and four rational 
numbers ai, where 240ai are nonnegative integer numbers such that 

ait'l + a2t2 + ^3^3 + Oiitl = /3fe, fc = 0, . . . , 7. 

Note that we have eight variables and eight equations. In fact, this problem 
has the unique solution: 



u 


-1 


-1/2 





1/2 


ai 


1 


56 


126 


56 



Geometrically we have the following picture: A„i(— 1, 1/2) n Am(y) = for 
y < 1/239, and for y = 1/239 the intersection of these convex sets consists of 
one point. In fact, this point is the vertex of the cone Am(j/). 

Step 3 : Feasible solution for a maximal code. 

It is well known fact that {vi} are vectors in R" if and only if its Gram 
matrix is positive semidefinite with rank that does not exceed n. 

Let C = . . . ,V24o} be vectors in and let G — (tij), tij = {vi,Vj), 
be the Gram matrix of this set of vectors. We have ones on the diagonal of G 
and other entries are —1,0, ±1/2. Let us denote by Qr{G) the total number of 
entries r in G. Now we have the following problem: 

To find a symmetric 240 x 240 positive semidefinite matrix G with ones on the 
diagonal and with other entries —1,0, ±1/2 such that rank(G) < 8 and 

Q-i{G) - 240, Qo(G) = 30240, Q_i/2(G) = Qi/2(G) = 13440. 

Bannai and Sloane [4, Theorem 7] show that this problem has the unique 
solution. Namely, C is isometric to the set of minimal vectors in Eg. 

Note that if we apply this method for 24 dimensions, as well as for other 
spherical codes that were considered in 4;, then we also obtain optimal codes. 
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2. Let us consider this approach for kissing arrangements in four dimensions. 
If we apply the SDPq for m = 3, 4, 5, 6, 7, then we obtain: 

c*(3) c*(4) = 26, c*(5) = c*(6) = c*(7) w 25.5584. 

From this it follows that the cardinality c of 7r/3-spherical code C is at most 25. 
Let 

nrniy) A™(-1, 1/2) n A„(y). 

As above, flrn{y* {m))) consists of one point. However, in four dimensions this 
point is not the vertex of Am{y)- 

3. If codes have some additional symmetries then the problem of bounding of 
codes sometimes can be solved easier than in general case. 

Let us consider spherical antipodal codes. We say that a spherical code C is 
antipodal ii v G C implies —vGC. Note that the kissing number problem for 
antipodal codes in dimensions n < 8 was solved by Boyvalenkov 

It is easy to see that for an antipodal code C we have 

1 

1 — c 

So in this case the independent variables are y and X2,X4, . . . , X2m~'2- 

If we apply the SDPq in four dimensions with to = 4, then we obtain 
c* = 24, 5/23, xl = 2/23, xl = 5/92, and x*^^^^ = -1/23. That 

yields (by the same method as in Step 2): 



u 


-1 


-1/2 





1/2 


Oil 


1 


8 


6 


8 



In fact, this implies the unique solution for the maximal code C: C is isometric 
to the set of minimal vectors in D4. 

In the same way maximal antipodal codes in dimensions 6,7 can be found. 
For n = 5 the method gives c* = 42. However, we expect that c = 40. It is 
possible using the fact that cai are integer numbers to prove that c < 40 (see 
details in [7]). 

6 An extension of Delsarte's method 

In this section we show how some recent extensions of Delsarte's method can 
be reformulated as a semidefinite programming problem. 

Let C = {zi} be a finite subset of a compact 2-point-homogeneous space M, 
and let 

Sk{C\[u,v)) -.^ ^5- 
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The constraints (3), (4) are relatively strong and sometimes give tight upper 
bounds for codes by Delsarte's linear programming method. However, there are 
some additional constraints for S'-codes which are not follow from (3), (4). 

For instance, Odlyzko and Sloane [H] in dimension 17 made use of the fol- 
lowing additional inequalities for a spherical 7r/3-code (a kissing arrangement): 

/3 

^at<l, ^at<2, h^-^, <2 = -^VS, 

t<ti t<t2 

i.e. 

so{C, [-1, h)) < \Cl so(C, [-1, h)) < 2\C\. 

These inequalities easily can be derived from the following fact (see [HE] or [25]): 
the cardinality of a spherical 9-code on an open spherical cap with angular radius 
ifk is at most k, where 

cosiy9fc = tk := v^cosS + (1 ^ cos9)/{k + 1). 

This fact for a spherical 9-code C yields 

so(C, hl,ife)) <k\C\. 

Actually, any upper bound on the cardinality of a spherical 9-coAe in a 
spherical cap or a spherical strip implies some inequality for so(C, [u,w]). In 
[5] an overview is given of the known bounds and several new upper bounds 
for codes in spherical caps and strips have been proposed. Agrell, Vardy, and 
Zeger [1] applied bounds on so(C, [u,w]) to improve some upper bounds for 
constant-weight codes. 

Recently, Pfender |22j found a new inequality for S2'- 

S2(C, [-1, -%/cos^)) < so(C, [-1, -%/cos^)) cos 61 -I- |C|(1 - cosS). 

Using this, he improved the upper bounds for the kissing numbers in dimensions 
9, 10, 16, 17, 25 and 26. Note that Krasikov and Litsyn [13] using a similar 
approach improved some upper bounds for binary codes. 

In the papers [THl dZl [T51 [TO] we found a few inequalities for some linear 
combinations of Sfe(C, [— l,to]) (0 < fc < 9), where C is a spherical 7r/3-code 
and <o < —1/2. In particular, that allowed us to prove that the kissing number 
in four dimensions is 24. This also yields a new solution of the Newton - Gregory 
(thirteen spheres) problem [17j . 

Let C be an S-code on M. As above, we assume that 5* = T n [a, 6] and 
$fc(i) is a polynomial of degree k. 

Let o = ui < 1*2 < . . . < tif < Uf-i-i = b and let for fc = 0, 1, 2 . . . 

xf' -.^ Sk{C,[ui,Ui+i)), xf^ -.^ Sk{C,[ui,b]). 

Then we have 

z^c + x^^'^ + ... + xf\ where z := |Cp, c := |C|. 
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In fact, the equality z = c? cannot be written as some positive semidcfinite 
constraints. However, the inequahty z > c? can be expressed in the form: 

r(c,.):=(^ ^)tO. (5) 

Using these notations (3) can be expressed as: 

k I 

c + ^Pfcd^xf^ >0, fc = l,2,... (6) 

d=0 i=l 

Since x^*^^ are power sums of distances Tq^r in the interval we can 

apply Theorem 3.2 to this interval: 

Hm{xf\xf\...,xf"'-^\[ui,Ui+r])^0, i = l,...,L (7) 

Note that all additional inequalities that were considered above can be ex- 
pressed in c, {xf^^}. For instance, let 

1 = 2, ui=a = —l, U2 = — VcosO, us = b = cos6. 

Then Pfender's inequality can be expressed in the form: 

x^i^ < cos 6* + (1 — cos 9) c. 

This inequality is linear in c, {x^*^'}, i.e. it is a positive semidefinite constraint. 

Suppose for a given £ and {ui} we have some additional positive semidefinite 
constraints for c and {a;-*^^}. Then we say that the set of these constraints is 
the set of Arf-constraints and write 

Ad(c,{a;f })tO. (8) 

Using (5) — (8) we introduce the following SDP problem: 
SDPA Problem: For given £, {ui}, Ad-constraints, and m choose real numbers 

fe=0,...,2m— 1 

X, > SO as to 

' ii=i,...,e 

maximize c 

subject to (5) — (8). 

Denote by c* (ni) the maximal value of c that can be attained in this problem. 
If C is an S'-code then c and xl certainly satisfy the constraints in the SDPA 
problem. Therefore c*(m) is an upper bound for j4(M, S). 

Theorem 6.1. For any natural number m we have 

A{M,S) < c*(m). 
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It is easy to see that c*(m + 1) < c*{m). Let us denote the mimmum of 
c*(m) by SDPA{M,'S>, [a,b],m,Ad). 

Theorem 6.2. SDPA{M, $, [a, b],m, Ad) < SDPo{M, $, [a, b],m). 

Proof. Note that the variable y in the SDPo problem can be expressed as 

c 

y = ■ 

z — c 

If we consider the constraints (6), (7) in the SDPA problem for the variables 
y,{xi} of the SDPq problem, then we obtain the constraints (3), (4). There- 
fore, the constraints in the SDPA problem imply the constraints in the SDPq 
problem. (Moreover, we have the additional constraint (8).) Then we have 

r* 



z — c 

2 



We also have (5): z> c which yields 

1 c* 



c* - 1 {c*f 

Thus 



> > V*. 



□ 

Clearly that the SDPA problem with empty Ad-constraints set has the same 
constraints as the SDPq problem. Then 

SDPA(M,<i>, [a,6],m,0) = SDPoiM,<i>, [a,b],m), 

i.e. the SDPA bound is equal to the Delsarte linear programming bound. 
Therefore, the SDPA bound with nonempty Ac?-constraints set can be con- 
sidered as an extension of Delsarte 's method. We see that for some cases this 
extension gives better bounds than the classical Delsarte method [13l [161 [El 
[I8jl9,22,- However, it is not easy to find (and to prove) new Ad-constraints. 



7 SDP bounds for codes 

In this section we explain how to use the p.d.k. property to obtain upper bounds 
on S'-codes in a subset 11 C M and to improve the SDPq and SDPA bounds. 
For any finite subsets X,Y in a 2-point-homogeneous space M we define by 

Fk{X,Y):= J2 '^k{rix,y)). 
It is easy to see that Theorem 2.1 implies the following fact: 
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Theorem 7.1. Let be finite subsets of M. Then for any integer 

k >0 the matrix {Fk[Xi, Xj)) is positive semidefinite. 

7-A. The SDPA bounds for codes on subsets. Let IT bo any subset of M. 
The largest size of an S-code in 11 will be denoted by ^(n, S). (As above, we 
assume that S = T r\[a,b].) 

Let Co = {qi, . . . , qr} be a fixed subset of M. Let — oo < ai < A < oo for 
i = 1, . . . ,r. We define IIi = ni({aj, by the following constraints: 

Hi = {a; e M : aj < t{x, qi) < i = l,...,r. 

Clearly, if aj = — oo, (3i — oc' for all i, then Hi = M. 

Suppose n C III Consider any 5-code C on 11. Let Xi = {qi}, i = 1, . . . ,r, 
and Xr+i = C. Then Theorem 6.1 yields 

Gk:={Fk{Xi,Xj))>Q, A; = 0,1,.... (9) 

Note that in the matrix Gk for 1 < i,j < r entries {Gk)ij = ^k{T{qi,qj)) 
are given. {Gk)r+i,r+i = Fk{C,C) (see (6)). Let 

qec 

Then y[°^ = y^°^ = . . . = y^P^ =c=\C\ and 

k 

{Gk)i,r+i = ^Pkdyf\ i < r. 

Thus, the matrix Gk linearly depends on the variables yj'^'' and a;^'^' . 
On the other hand, since ai < T{q, qi) < (3i, q G C, we have 

iJ„(2/f\2/«,...,2/f"-^\[a„/3,]))^0, ^ = l,...,r. (10) 
For a given m let us consider the following SDP problem: 

maximize c 

s.t. 

(7), (8), (9), (10). 

Denote by c*(m) = c*{m,{qi},{ai},{(ii}) the maximal value of c that 
can be attained in this problem, and let us denote the minimum of c*{m) by 
SDPACni,^, [a, b], Ad). Then we have the following theorem: 

Theorem 7.2. ^(n,^) < SDPA{Ili,^, [a,b],Ad) < c*(m) 
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7-B. The SDP bounds. Now using Theorem 7.1 wc improve the SDPA 
bound. Let Co = {qi, . . . , Qr} be a fixed subset of M and let 11 C M. Suppose 
a family {Ili}i=i,...,d is a partitioning of 11, 

n c Hi u . . . u Hd, Hi n n, = 0, i j. 

We also assume that 

C n^- := {g e M : aij < T{qi,q) < Pij, i = l,..., r}. 
For finite C C 11 denote by 

d—IlinC, Ci:=\Ci\. 



(k) 



^^r'^K?;), l<i,j<d, fc = 0,l,2,..., 

ueCi veCj 



•= E ^'^ l<^<^. ^<j<d, fc = 0,1,2, 

qeCj 

Recall that ^k{t) = J2Pkit'- Then 

fc 

Fk{Ci,Cj) = J2 Pu ^1? , i^fc (ft , Q ) = ^ pu y\f . 

i=0 1=0 



If we apply Theorem 7.1 with Xi = {qi\,i = 1, Xi+r = Ci,i = 

1,. . . ,d, then we obtain that the matrix Fk = {Fk{Xi,Xj) is positive semidef- 
inite. We see that entries of Ffe are linearly dependent on the variables x = 
jx-^^l and y = (Obviously, a;^-^^ = x\j\ So we can assume that x is a 

set of variables a;-^^ with i < j.) Therefore, Ffc = F^,(x, y). Thus for any m > 0: 
Ffe(x,y)^0, fc = 0,l,...,2m-l. (11) 
Note that all a;^°'' , y^j^ are nonnegative integer numbers and 

yS^=cj, xf=c,cj. (12) 
Since aij < T{q,qi) < pij, q G Cj, we have 

Hmivlf , vlj^ yg""'^ > Ki ,Pij])hO; i = l,...,r; j = l,...,d. (13) 

Let C be an S'-code on 11. All values {t{u, v)},u E Ci,v G Cj,u v belong 
to some interval [aij,bij] C [a,b]. It is clear that the bounds {aij,bij} depend 
only on 11^,11^, and [a,b]. 

By assumption for i ^ j we have 11^ n Ilj =0 which yields C, fl Cj = 0. 
Then t{u, v) S [a^ , bij] for all u G Ci,v G Cj. Therefore, 

Hm{x\f,x\]\...,xl^.'"-'\[aij,bij])hO, l<i<j<d. (14) 
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Note that among ordered pairs u,v G Cj wc have exactly q pairs 

{u,u) such that t{u,u) = tq ^ [dajbu]. All other pairs belong to this interval. 
Let Zk,i := x-f' - c, Tq , fc > 0, 2o,i := — Cj. Then we have 

Hmizo.i,Zi,i,...,Z2m-l.i,[o-ii^bii])yO, i=l,...,d. (15) 

In Section 6 we considered some additional constraints which don't follow 
from (3), (4). In the case when with (11), (13), (14), (15) we have some additional 
positive semidefinite constraints we write them in the following form: 

Ad{^,y,{ni},[a,b])yO. (16) 



Now we introduce the QOP bounds for codes on 11. 

QOP Problem: For given m, Co, d, {Hj}, and Ad-constraints choose nonnegative 
integer numbers Cj and real numbers x, y so as to 



maximize c = ci 



subject to 



(11) - (16). 



Actually, the QOP problem is not a semidefinite programming problem. In 



(12) we have (z) 



JO) 



iCj, and (m) Ci are nonnegative mte^er numbers. So {i) 



and (ii) are not positive semidefinite constraints. 

Let us consider a and xf^^ as independent real variables with the following 
positive semidefinite constraints: 



(0) 



Cj > 0, 



/ 1 

Cl 



Cl 

/ill 



V Cd hid 



Cd \ 

hid 
hdd J 



hij 



„(0) 



(17) 



If we substitute (17) with real Cj for (12) in the QOP problem, then we get a 
semidefinite programming problem. Let us denote this problem by SDP. 

Denote by c*{m,Co,d,{Ili},Ad) the maximal value of the sum c that can 
be attained in the SDP problem. Then 



Theorem 7.3. A{n,S) < c*{m,Co,d,{Ui},Ad). 



7-C. New upper bounds for one-sided kissing numbers. Our first com- 
puter experiments show that this method can significantly to improve some 
upper bounds for codes. To show that we consider the one-sided kissing prob- 
lem. 

Let H he a, closed half-space of R". Suppose 5 is a unit sphere in H that 
touches the supporting hyperplane of H. The one-sided kissing number B{n) is 
the maximal number of unit nonoverlapping spheres in H that can touch S. 



17 



If nonoverlapping unit spheres kiss (touch) the unit sphere S in H C R", 
then the set of kissing points is an arrangement on the closed hemisphere S+ 
of S such that the (EucUdean) distance between any two points is at least 1. 
So the one-sided kissing number problem can be stated in another way: How 
many points can be placed on the surface of S+ so that the angular separation 
between any two points is at least 7r/3? In other words, B{n) is the maximal 
cardinality of a 7r/3-code on the hemisphere S+. 

Currently B{n) known only for n < 4. Clearly, B{2) — 4. It is not hard to 
prove that B{3) — 9. Recently, using some extensions of Delsarte's method we 
proved that B{4) = 18 (see [12] for a proof and references). 

In [51 Section 7] we derive estimates of B{n) for n = 5, 6, 7, 8. Actually, these 
bounds are based on the LP bounds for kissing numbers and the LP bounds on 
^o-codes, where cos^o — l/\/3- So we call these upper bounds the LP bounds. 

We applied the SDP with [a,b] = [—1,1/2], the empty Ad - constraints set, 
m = 6, r = 1 (i.e. Co consists of the one point eo), d = 3, and 

ni-{peS"-i:dist(p,eo)<V4}, 

Ha = {p e S"-i : 7r/4 < dist(p, eo) < 2tt/5}, 
Ha = {p e S""^ : 2n/5 < dist(p, eo) < tt/2}. 
We give the SDP (as well as LP) bounds on B{n) in the following table: 



n 


3 


4 


5 


6 


7 


8 


9 


Lower bound known 


9 


18 


32 


51 


93 


183 




LP bound 


9 


20 


39 


75 


135 


238 


378 


SDP bound 


9 


19 


35 


64 


110 


186 


309 



Note that for n = 8 the SDP bound is close to the known lower bound: 
B{8) > 183. So using some new Ad - constraints or other families {11^} we have 
a good chance to solve the one-sided kissing problem in eight dimensions. 
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